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Abstract 

A semiclassical Thomas- Fermi method, including a Weizsacker gradient term, 
is implemented to describe ground states of two dimensional nanostructures 
of arbitrary shape. Time dependent density oscillations are addressed in the 
same spirit using the corresponding semiclassical time-dependent equations. 
The validity of the approximations is tested, both for ground state and density 
oscillations, comparing with the available microscopic Kohn-Sham solutions. 
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I. INTRODUCTION 

Two-dimensional quantum dots, and semiconductor nanostructures in general, are exam- 
ples of artificial systems in which, ultimately, we may seek to produce and control electronic 
quantum properties. Sometimes they are referred to as artificial atoms and in fact, it has 
been proved that the electronic structure of circular quantum dots resembles in some as- 
pects the shell structure found in atoms and nuclei. Particularly relevant in this sense are 
the Coulomb blockade measurements that provide direct acces to the electronic energy levels 
by adding electrons one by one to the quantum dot [^ . The magic numbers measured in ver- 
tical quantum dots [Q, corresponding to maxima in the addition energies, have also shown 
clear evidences of shell structure. Absorption measurements in the far-infrared region j^,^ 
as well as experiments using sophisticated light scattering techniques [^^ have probed both 
charge density (CDE) and spin density (SDE) excitations and have shown that parabolic 
confinement is closely satisfied for small size dots. 

Interest in the quantum dot community has recently focussed on the properties of de- 
formed nanostructures. For instance, elliptic quantum dots have been investigated in Refs. 
|7Hr3| . Particularly, in Refs. [lllJl^ some of us have addressed the collective oscillations of 



deformed dots where, in addition to density and spin modes, it has been predicted the exis- 
tence of orbital current modes at low energies. This interest in deformed dots is motivated 
by the advances in nanofabrication techniques, that allow to produce quantum dots of many 
different shapes. 



Microscopic theoretical approaches, like Hartree [Q , Hartree-Fock ||T5| and density func- 



tional HlBHlSJ, let alone exact diagonalizations |jT9|,^, get very demanding computationally 



for increasing number of electrons, especially in the symmetry unrestricted case. It is thus 
interesting to develop a semiclassical approach, dealing only with the total density, for which 
the computational effort is not much dependent on the number of electrons. This is our pur- 
pose in this paper. Semiclassical Thomas-Fermi models have been already used in the field. 



for instance in Refs. [pl|-p^. In particular, Ref. [^] provides a quite rigorous presentation 
of the theory. However, previous works deal only with the semiclassical ground state and 
mostly for circular symmetry. Here we emphasize the application to deformed structures 
and concentrate on the oscillation modes. We do not include in this paper magnetic field 
since this requires a non trivial generalization of the theory that we plan to develop in a 
separate contribution. 

The structure of the paper is as follows. In Sec. II we discuss the ground state in our 
semiclassical approach and compare with available microscopic Kohn-Sham (KS) results. 
Section HI introduces the formalism for the time-dependent oscillations. Results for the 
oscillation modes are discussed in Sec. IV. Finally the conclusions are drawn in Sec. V. 

II. GROUND STATES 

A. Definition of energy functional 

Using Density Functional Theory in a local approximation we write the total energy in 
terms of the electronic density p(r) as E[p\ = J drS[p], with the energy density separating 
as 



8[p] = t[p] + -VH{r)p + Sxcip) + Vext{r)p ■ (1) 

The different terms are the kinetic energy density t[p], the Hartree potential t'_f/(r) = 
/ dr' f_J. , the exchange-correlation energy density Sxcip) ^^d the external confining poten- 
tial fe3;t(r)- The Thomas- Fermi approximation in two dimensions (2D) yields the following 
kinetic energy density |^3[ 

-Hp1-£(V + ^vv). (2) 

By analogy with the Weizacker term for the 3D kinetic energy functional, we have added a 
gradient term that gives the exact kinetic energy for a single electron 

rwip] = ^A-— , (3) 

with A = 1/4. Our total Thomas-Fermi- Weizsacker (TFW) kinetic functional is then r = 
Ttf + Tw- It is worth to point out that the first no n- vanishing gradient correction in 2D is not 
known from a rigorous semiclassical expansion and therefore we have introduced empirically 
the Weizsacker-like term. We will show later that the results are not sensitive to the precise 
value of the coefficient A. For the sake of comparison, we recall here that the KS method 
provides the exact kinetic energy by means of a set of single particle orbitals {v^i} as 

r(r) = ^ Y. \^V^i^)\" ■ (4) 



2m 



l,OCC. 



For the exchange-correlation energy £xc we have used the same functional of Ref. |jT^, based 



on the Tanatar and Ceperley calculation for the uniform gas EH]. Specifically, in modified 



atomic units 27 , it is written as 



r^ / \ 4 /2 o/n 1 1 + aix , . 

^xc{p) = -^\-p' + o^'^Ptt 1 ^T z ' 5 

6\t: 2 1 + aix + a2a;^ + a^X'^ 



where x = (np) ^1^ and the a^ coefficients are given in Ref. |2^ . Notice that in this formalism 



we always assume perfect spin degeneracy in the ground state in order to obtain a single 
chemical potential for both spin components (see below). 

The ground state density is determined from the minimization (Euler-Lagrange) equa- 
tion, with the contraint of particle number conservation 



^ E\p\ - p f drp) =0 . (6) 



6p 



Specifically, this reads 



2m 2m, p 



+ ( Vext + vh + —^ + '^^l^P ] P = PP ■ (7) 



By making the transformation ip = y/p this equation may be written in the famihar form of 
a Schrodinger hke equation 

- 4A-— V'^ + F^ = p^ , (8) 

where we have defined the semiclassical effective potential V = f ext + ^H + ^^^ + Svr^p. The 
chemical potential /x plays the role of the eigenvalue in Eq. (H). Written in this way, we may 
now use the algorithms developed to solve the KS equation in arbitrary 2D confinements 
| rT| . Our method is based on a discretization of the xy plane in a grid of uniformly spaced 



points. Then, the total number of grid points, not the number of electronic orbitals as in KS 
theory, determines the computational cost of the problem. Typically, we use grids ranging 
from 50^ to 100^ points. 

B. Results and comparison with Kohn-Sham 

In this subsection we show results for the semiclassical ground states in different confin- 
ing potentials, focussing especially on the comparison with the corresponding microscopic 
KS results in order to prove the validity of the approximations. Since the microscopic calcu- 
lation for medium and large sizes is computationally feasible only when circular symmetry is 
imposed to the system, we will begin by considering the cases of circular parabolas and disks 
of jellium. In these cases we compare with the KS radial solution |]T^. Then we analyze 
one case with deformation, namely the deformed parabola with A^ = 20, and compare with 
a symmetry unrestricted microscopic calculation |]Tl] . 



1. Circular parabolas 



The confining potential in this case is (we will use modified atomic units [^] for the rest 
of the paper) 

Ve,t{r) = Ao + -ujy. (9) 

The constants Aq and ujq are usually parameterized to reproduce the potential and its 
curvature at the origin of a jellium disk (see below) with Np positive charges as Aq = 
—2J Np/rs and uJq = l/{JNpr^), with r^ the radius per unit charge. Figure 1 shows the 
densities obtained taking r^ = l.Slag. Panels (a) and (b) prove that the semiclassical density 
closely adjusts to the microscopic one for varying electron number in a fixed parabola and 
for a varying parabola curvature at fixed electron number, respectively. As is well known 
from atomic physics, the TFW density averages the microscopic shell oscillations in the inner 
region. Panel (c) shows that the coefficient A controls the surface width of the semiclassical 
density. The value A = 1/4 provides a reasonable approximation to the KS density tails in 
all cases. However, it is clear that if this coefficient is allowed to vary a better fit of the 
densities may be obtained. We have not followed this procedure since the fitted value of A 
would be different for each type of confining potential (see next subsections). 



In Fig. 2 we compare total energies per electron in the two approaches. The fit of 
energies is amazingly good; the deviation of the TFW energies with respect to KS being 
less than 0.5% for all cases shown in Fig. 2. Although the high precission of the TFW 
model for this quantity may seem a bit surprising, we recall that the method looks for 
the variational minimum of total energy. Therefore, this is a priori the best value of the 
model. We have made the comparison for magic electron numbers, corresponding to closed 
shell configurations. This is the reason why shell oscillations with size are not visible in 
the microscopic results of Fig. 2(a). The total energies are almost independent of A for 
reasonable values of this coefficient. The results for iV = 42 and A = 0.25, 0.5 and 1.5 are 
E/N = -4.033, -4.025 and -4.014 H*, respectively. 

2. Jellium disks 

Another external potential with circular symmetry is that of a uniformly charged disk. 
Defining the positive charge density (jellium density) in terms of the r^, parameter pj = 
l/(7rr^) and for a disk with radius R, the potential is 



Vext{r) 



4 R 

— ^E{r/R) if r <R 



4 r 



E{Rlr)-[\-(?-^\^K{Rlr) 



(10) 
if r>R 



where E and K are the elliptic integrals. Assuming a uniformly charged disk the number of 
positive charges A^j, is related to the disk radius by i? = rgJNp. 

Contrary to the parabola, the jellium potential forces the electronic density to saturate 
inside the dot. This is obviously due to the charge screening effect, that energetically favors 
the cancellation of charges. Figure 3 shows the electronic densities of neutral disks with 
Tg = 1.51aQ as a function of size. The agreement between both models is rather good. As 
for the parabolas, the TFW densities average the oscillations of the KS ones but now they 
rapidly saturate inside the dot. Close to the edge, the TFW densities present an small 
oscillation. This is similar to the Friedel-type oscillations found in metal surfaces and is 
enhanced by potentials that abruptly vanish at the edge, as the jellium one. Also shown in 
Fig. 3 (panel b) are the energies per electron of the neutral disks, that are also very well 
reproduced by the semiclassical model. 

From Fig. 3 we see that the TFW density tails slightly overestimate the KS ones, with 
the used value of A. This implies that the number of electrons that are spilling out of the 
jellium disk will be slightly enhanced in the semiclassical method. This spill out mechanism 
is known from cluster physics to be of fundamental importance for a proper description of 
the optical properties [^ . The fact that our TFW method correctly includes spill out gives 
us some confidence in its use for the description of time dependent density oscillation (see 
Sec. IV). 

3. Deformed parabola for N=20 
The third case we have considered is an elliptical dot, confined by an anisotropic parabola. 



VeAr) = ^-u',j^^^{x' + f3V). (11) 

The parameter (3 gives the ratio of parabola coefficients in y and x directions, i.e., writing 
the external potential as Vexi(r) = |(<^x^^ + ^lu"^)^ we have [3 = Uy/oj^- At the same time 
the centroid {ujx + ujy)/2 is kept fixed at the value cuq, defined as for the circular parabola 

^l = l/(-?v^)- 

Since the KS problem is much more involved than for the preceding circular potentials, 

we restrict here the comparison to the A^ = 20 electrons case, with Np = 20, and consider 

four deformations (3 = 0.875, 0.75, 0.625 and 0.5. Figure 4 shows the densities corresponding 

to these four deformations as well as the circularly symmetric case {f3 = 1) for completeness. 

It is seen that, quite nicely, the elliptic contour lines of the TFW density follow in average the 

equidensity regions of the KS result. For /? = 1 and /3 = 0.875 the structure of equidensity 

regions of both models is very similar. For lower values of (3 the microscopic model yields an 

incipient electron localization that the semiclassical model is obviously not able to reproduce; 

it gives nevertheless the correct average value. We conclude that the semiclassical model 

reproduces the average density distributions also in non circular systems. Total energies, 

given in Tab. I, are also in excellent agreement with the KS ones. 



III. TIME DEPENDENT EQUATIONS 

To derive the time dependent equations in the semiclassical approximation we will follow 
the fluid dynamical approach. This method has been used in nuclear physics [^ and, more 



recently, also for electronic oscillations in atomic clusters [^. The variational derivation of 
these equations can be found in the book by Ring and Schuck [^. Here we will just point 
out the essential ingredients and particular details for the application to our case. 

The essential assumption is that all the single particle orbitals evolve in time with a 
common complex phase as 

y,,(r,t) = y,f(r,t)e-(^'*), (12) 

where both </?,• (r, t) and s(r, t) are assumed to be real functions. With (0) the kinetic 
energy separates in two contributions 

T = Tintr + -I drpu^ , (13) 

an intrinsic one Tj„fr = ^J2i \^^i P and another associated to the common velocity field 
u = Vs. Generalizing this separation we may write for the total energy (expectation value 
of the Hamiltonian H) 

n={H) = E,ntr + \jdrpv? . (14) 

In our density functional approach the intrinsic energy is in fact given by the energy func- 
tional Eintr = E[p\. Noticing now that p and s are conjugated canonical variables, Hamilton's 
equations are [^ p = ^, —s = ^. Specifically, these read 
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p = -V(pVs) 

-.4(V.)= + ^. (15) 

The two Eqs. (|T^) yield the time evolution of the semiclassical variables p and s, and thus 
they are our required input to model a semiclassical dynamics in quantum dots. As happens 
for the ground state equation, we may still transform Eqs. (|l^) to a form similar to the 
microscopic equation. In fact, defining the time dependent complex function ip = y^e'* 
both Eqs. ([ISD reduce to 

i.e., an equation identical to the time-dependent KS equation if we identify the contribution 
within brackets in the right hand side with the potential term. 

In the preceding discussion we have assumed, for simplicity, that both spin densities 
are oscillating in phase, and thus identically to the total density. The formalism, however, 
can also account for the situation in which both spin densities p^ (77 =1, i) oscillate out 
of phase, as happens for instance in spin modes. In this case we just need to assume the 
semiclassical approximation for each spin fluid, with the Coulombic coupling terms between 
them and generalize the exchange-correlation contribution to the locally polarized case using 
the exchange interpolation formula [|26|. Defining ipj^ = ^p e'^^^ the two equations are 

In this paper we apply the time dependent equations to obtain the linear response fre- 
quencies corresponding to dipole charge density (CDE) and spin density (SDE) excitations 
of general 2D nanostructures. These are normal modes of oscillation and the technique to 
obtain them is simply to apply an small instant perturbation on the ground state and take 
this as initial condition for the time simulation. For the CDE the initial perturbation is 
simply a rigid translation of the total electronic density, given by operator T(a), with a 
the vector displacement. For the SDE we apply opposite translations for spin up and down 
densities, 7^(a^), with a^ = — a^. After this, we keep track of the time evolution of the dipole 
moments (i^(t) = {r)n ■ e, where e is a unitary vector in the direction of the displacement, 
and 

(r), = ^/rfrp,(r,t)r, (18) 

with N^, Ni the number of electrons with spin up and down, respectively. A frequency 
analysis of the total dipole moment d^ + di for the CDE, and of the spin dipole moment 
(i| — di for the SDE, finally yields the oscillation frequencies 



IV. RESULTS OF DENSITY OSCILLATIONS 

We present in this section the spectra obtained within the semiclassical formalism. We 
consider density and spin dipole oscillations and, as in Sec. I, we emphasize the comparison 



with the corresponding KS results. One case is taken as representative for the three types 
of confining potentials for which we have already discussed the ground state. 

4. Circular parabola with N = 56 

Figure 5 shows the CDE and SDE for the parabola with A^ = 56 and Np = 40. The 
results have been plotted in a logarithmic arbitrary vertical scale (each tick marks an order 
of magnitude). First, we notice that in the CDE the semiclassical spectrum (and also the KS 
one) yield a single frequency that coincides with the parameter uq of the confining potential. 
We thus conclude that time-dependent TFW satisfies well the generalized Kohn's theorem 
p!9| . This theorem states that the exact dipole CDE for parabolic confinement is a single 



peak at the parabola frequency luq. The more intense SDE's lie at lower energy because 



the residual interaction is attractive in this channel [0. The semiclassical model does not 
include the coupling to particle-hole excitations, that leads to an important fragmentation of 
the KS spectrum. Nevertheless, it reproduces some fragmentation of the collective strength, 
contained in two dominant peaks at ~ 0.05 H* and ~ 0.19 H*, that nicely correspond to 
very intense KS excitations. 

We remark that the microscopic results in Fig. 5 (and also those of Fig. 6 below) have 
been obtained by using the perturbative random-phase approximation (RPA) in a particle- 
hole basis fl^ while the TFW results were obtained, as explained in Sec. IV, from a frequency 
analysis of the real time signal. The RPA calculation provides a very high frequency reso- 
lution, even for low intensity peaks. On the contrary, the analysis of the time signal is not 
able in some cases to discriminate the low intensity peaks because of the limitations of a 
discrete time sampling and a finite total time window. 

5. Circular jellium with N = 58 

This is shown in Fig. 6. In this case generalized Kohn's theorem does not hold since the 
external potential is not of parabolic type. As a consequence the dipole charge oscillation 
couples to the relative motion and the spectrum is generally fragmented, with peak energies 



depending on the number of electrons |31[ . The disk potential behaves quadratically close to 
the disk center but it deviates for points closer to the edge. Fig. 6 shows the Uq value for the 
quadratic behaviour at the origin. The semiclassical model reproduces the blue shift from 
ujq although it slightly underestimates its quantitative value. This may be traced back to 
the slight difference in spill out mentioned in Sec. II. B. 2. The TFW model yields a greater 
spill out and thus produces a softer oscillation mode. The absence of an exact Kohn mode 
manifests with an important fragmentation, partly reproduced in the TFW model. In the 
SDE of Fig. 6 the situation is similar to the circular parabola. The semiclassical spectrum 
reproduces the dominant collective peaks, and has less fragmentation than KS. 

6. Deformed parabola 

Figure 7 shows the results for an anisotropic parabola with A^ = 20, A^p = 20 and different 
deformations. The two upper panels correspond to the CDE for j3 = 0.5 and 0.75. They 
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prove that also in deformed parabolas the generalized Kohn's theorem is well satisfied by 
time dependent TFW. In this case the parabolas in x and y directions are different and, 
in fact, the oscillation in these two directions is at frequency uj^ and Uy, respectively. For 
the two other deformations not shown in the figure (/3 = 0.625, 0.875) generalized Kohn's 
theorem is equally satisfied. Therefore, as for TDLDA, TDTFW fulfills the exact property 
for dipole charge oscillations saying that the two center of mass coordinates X = Y^iXi and 
Y = J2iyi oscillate with the frequencies Ux and ujy, respectively 



The four lower panels of Fig. 7 show the SDK spectra for the four deformations considered. 
The x-y splitting in the spin channel is nicely reproduced, as a function of deformation within 
the TFW model. Fragmentation is present in both models, although TFW is overestimating 
the collective strength at high energy for /? = 0.75 and (3 = 0.875. We attribute this to the 
rather small electron number of this dot {N = 20), for which TFW is surely less accurate 
than for large sizes. We remark again that the lower energy (and more intense) peak, as 
well as the x-y splitting are correctly reproduced. 

To finish this section and in order to emphasize the power of the semiclassical method 
we reproduce in Fig. 8 a case that is beyond the present capability of microscopic KS 
calculations. It corresponds to 72 electrons in a deformed parabola, with Np = 72 and 
Tg = 1.51aQ. Since the CDE corresponds only to u^ and Uy, according to Kohn's theorem, 
we display the result of the SDK for (3 = 0.75, as well as the circular one for comparison. 

V. CONCLUSIONS 

In this paper we have discussed a semiclassical approach to the ground state and density 
oscillations of 2D nanostructures. The method has been implemented for the general case 
in which no spatial symmetry is required. The validity of both ground state and dipole 
oscillation descriptions has been checked by systematically comparing the semiclassical TFW 
results with the corresponding KS ones. We have shown that the TFW densities closely 
follow the KS ones, averaging the shell oscillations, for different types of external confinings: 
circular and deformed parabolas and jellium disks for which the density saturates. Besides, 
the TFW model reproduces the KS energies with great accuracy. The dependence of the 
energy with the value of the Weizsacker coefficient is very small. This coefficient controls 
the density tail and electronic spill out for jellium disks. We have shown that the value 
A = 1/4 provides a good overall fit. 

Dipole charge density and spin density oscillations have been analyzed by using the 
time dependent semiclassical equations. In circular parabolic confinements we have shown 
that the semiclassical spectra satisfy well the generalized Kohn's theorem for charge density 
excitations, while in neutral jellium disks it yields a blue shift similar to the KS one. In 
elliptic dots the generalized Kohn's theorem is satisfied as well. For spin density excitations 
the TFW model is able to reproduce the dominant peaks of the spectrum and the splitting 
associated with deformation in elliptic dots. 

In general, we have shown that the semiclassical Thomas-Fermi- Weizsacker model pro- 
vides a reliable tool to obtain accurate approximations to the ground state and linear os- 
cillations of quantum dots. This opens the possibility to use it in order to explore a great 



variety of confining potentials witli different geometries and for large sizes. 
This work was performed under Grant No. PB95-0492 from CICYT, Spain. 
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FIGURES 

FIG. 1. Radial electronic densities for circularly symmetric parabolas in the TFW model 
(dashed) in comparison with the corresponding KS densities (solid). Panel (a) corresponds to 
a fixed parabola with ujq = 0.21 H*, that with the parameterization used in the text is given by 
Np = 40 and r^ = l.Slag, with different number of electrons. Ordered vertically at r = are shown 
A^ = 2, 6, 12, 20, 30, 42, 56, 72 and 90, respectively. Panel (b) corresponds to a fixed number 
of electrons A^ = 42 and varying parabolas with the same r^ and Np = 20, 30, 40, 50 and 60 (also 
ordered vertically at r = 0) . Panel (c) shows the density for a fixed dot, iV = 40 and Np = 42 
within KS and for varying Weizsacker coefficient A = 0.25, 0.75 and 1.5 (in vertical order around 
r = 14). 

FIG. 2. Energy per electron, as a function of electron number A^, in the TFW model (smooth 
line) and in the KS model (symbols). Panels (a) and (b) contain the results for the corresponding 
panels in Fig. 1. 

FIG. 3. Results for neutral jellium disks with r^ = 1.51aQ. Panel (a) shows the radial densities 
within the TFW (dashed) and KS (solid). With increasing edge radius the results correspond to 
N = 12, 20, 46, 58 and 72. Panel (b) shows the corresponding energies for KS (symbols) and TFW 
(smooth line). 

FIG. 4. Comparison of KS and TFW densities for anisotropic parabolas with N = 20, Np = 20 
and Vg = 1.51ao for different deformation parameters. The elliptic contour lines correspond to the 
semiclassical densities while the coloured background regions show the KS results. The scale bar 
indicates the density associated with each colour as well as its value for the semiclassical contour 
lines, which correspond approximately to the colour transitions. 

FIG. 5. Dipole excitation spectra, in logarithmic arbitrary units, for charge density and spin 
density excitations in a parabolic quantum dot. In each panel, the upper curve corresponds to the 
microscopic RPA calculation and the lower one to the present semiclassical approximation. In the 
upper panel the triangle shows the value of the parabola coefficient cjq. 

FIG. 6. Same as Fig. 4 for a jellium disk with r^ = 1.51aQ. In this case the triangle corresponds 
to the parabola that fits the disk potential at the origin. 

FIG. 7. Dipole spectra for the parabolic quantum dots of Fig. 4. Solid and dashed lines show 
the results for oscillations in x and y directions, respectively. For the CDE (two upper panels) we 
plot the TFW result and the position of the tOx and tOy parameters while for the SDE (four lower 
panels) we compare with the corresponding KS curves (displaced vertically above the TFW ones). 
As before, vertical scales are in logarithmic arbitrary units. 

FIG. 8. Same as Fig. 7 for a parabolic dot with N = 72 and Np = 72. The microscopic KS 
results are not available for /3 = 0.75. 
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TABLES 

TABLE L Energies per electron (in H*) within the TFW and KS models. The results corre- 
spond to a parabola with N = 20, Np = 20, Vg = 1.51ag and different deformations. See text. 




/3 
1 0.875 0.75 0.625 0.5 


KS 
TFW 


-3.026 -3.032 -3.059 -3.118 -3.223 
-3.025 -3.033 -3.061 -3.119 -3.225 
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